Magnetic properties of frustrated spin ladder 



O 

o 
o 

(N 
00 



D 

■4— » 
•4-* 



c 

O 

o 



> 



(N 
O 
O 

o 

-I— > 



T3 

O 
O 



X 



Toru Sakai and Nobuhisa Okazaki 
Faculty of Science, Himeji Institute of Technology, Ako, Hyogo 678-1297, Japan 

(November 1999) 

The magnetic properties of the antiferromagnetic spin ladder with the next-nearest neighbor 
interaction, particularly under external field, are investigated by the exact diagonalization of the 
finite clusters and size scaling techniques. It is found that there exist two phases, the rung-dimer 
and rung-triplet phases, not only in the nonmagnetic ground state but also magnetized one, where 
the phase boundary has a small magnetization dependence. Only in the former phase, the magneti- 
zation curve is revealed to have a possible plateau at half the saturation moment, with a sufficient 
frustration. 

PACS numbers: 75.10.Jm, 75.30. Kz, 75.50.Ee, 75.60.Ej 



I. INTRODUCTION 

The frustration of the antiferromagnetic exchange in- 
teraction brings about many interesting phenomena in 
the quantum spin systems, because it generally enhances 
the quantum fluctuation. It would be valuable to con- 
sider the effect of the frustration on the spin ladder, 
like the materials SrCu 2 3 (Ref. §), CuafQjH^Na^CLi 
(Refs. |,§) and La 6 Ca 8 Cu240 4 i (Ref. {§). They are 
strongly quantized and have the spin gap. When the 
next-nearest-neighbor (NNN) exchange interaction ap- 
pears, the frustration takes place in the system. In the 
classical limit it is easily shown that the system has two 
different ordered phases depending on the strength of 
the NNN exchange and the phase boundary does not 
change even under external magnetic field. In the quan- 
tum system, however, some modifications should exist in 
the ground state phase diagram, because the spin ladder 
has no long range order even at T = 0. In this paper, we 
investigate the frustrated spin ladder by the exact diago- 
nalization of the finite clusters to determine the magnetic 
phase diagram, even under external field. In addition we 
consider the possibility of the magnetization plateau || , 
which is predicted by a strong coupling approach. M 



II. MODEL AND NUMERICAL METHOD 

The S — 1/2 spin ladder with NNN coupling is de- 
scribed by the Hamiltonian 

L 

H. = Ji / ]{8i t i ■ Si^+i + S2,j • S2,;+i) 

i 
L 

+ Jj_ X)( S M • S 2,i) 

i 
L 

+ J2 / j(Si,j • S2,i+l + S2,i • Si^+i), (1) 

i 

where J\ , J2 and J± are the coupling constants of the leg, 
NNN (diagonal) and rung exchange interactions, respec- 



tively. We put Jj_=l in the following. Using the Lanc- 
zos algorithm we numerically solved the ground state of 
the finite clusters. We also calculated the lowest en- 
ergy of H for J2i( S i,i + s 2,i) = M > which denotes 
E(M). Using E(M), we investigate the magnetic state 
with m = M/L under the external field described by 
Wz = -^Ef(^, i + 5| ii ). 



III. TWO MAGNETIC PHASES 



Consider the nonmagnetic ground state at first. In the 
classical limit the system has two different ordered phases 
divided by the first-order phase boundary J 2 — J±/2(= 
1/2) shown as a dashed line in Fig. 1. Ji = 1/2 is also 
the boundary because the phase diagram should be sym- 
metric under the exchange of J\ and J2 (the reflection 
with respect to the dot-dashed line in Fig. 1). 
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FIG. 1. Phase diagram of the frustrated spin ladder 
with magnetization m for L = 12. Every phase boundary 
is first-order within this region. The dashed line is the classi- 
cal limit. The dot-dashed line is the symmetric line. 
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In the quantum S = 1/2 system we should distinguish 
the two phases based on the dimer picture; the dimers 
along the rung and the diagonal, respectively. The former 
is realized for J2 <C J±/2, while the latter for J2 3> Jj_/2. 
In the latter phase each two spins coupled by the rung 
are expected to behave like an effective S = 1 (triplet) 
object. Thus we call the two phases 'rung-dimer' and 
'rung-triplet', respectively. The phase boundary is easily 
detected as a level crossing point in the ground state even 
in small finite clusters. Since the boundary is almost in- 
dependent of L, we show only the result of L = 12 as 
circles in Fig. 1. Our study of the spin correlation func- 
tion along the rung also supported the above argument 
and suggested that the boundary is first-order. The re- 
sults are completely consistent with the recent analysis 
by the density matrix renormalization group. [Q (It also 
indicated the crossover of the phase boundary from first- 
order to second-order ones for J 2 < 0.287Ji, but we don't 
consider such a parameter region in this paper.) 
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FIG. 2. Canted Neel orders of the classical system under 
external field H (a) for J 2 < 1/2 and (b) for J 2 > 1/2. 

Even in the magnetic state under external field the two 
phases still can be identified by the different canted Neel 
orders shown in the Figs. 2(a) and (b), respectively, in 
the classical system. The phase boundary is the same as 
the nonmagnetic ground state. The quantum system is 
gapless for < m < 1 and it might be difficult to distin- 
guish the two phases by the dimer picture. In this case 
the classical picture is useful because the gapless phase 
is characterized by the power-law decay of the dominant 
spin correlation function corresponding to the classical 
order. Thus the quantum system should also have two 
phases like the classical limit. The same analysis as the 
nonmagnetic state indicated the first-order boundary for 
finite m. We show the boundaries for m=l/6, 1/3, 1/2, 
2/3 and 5/6 (£=12) in Fig. 1. They exhibit a small 
m dependence, although it is not so large that a field- 
induced transition between the two phases is expected 
to occur in any realistic situations. As the magnetiza- 
tion increases, the boundary tends to approach to the 



classical limit for most magnetizations. For m = 1/2, 
however, the boundary exhibits a quite different behav- 
ior and it is close to the nonmagnetic one. It implies that 
the quantum fluctuation is enhanced by the frustration 
particularly at m = 1/2. Thus we consider the possibil- 
ity of another spin gap induced by external field, that 
is observed as a plateau in the magnetization curve at 
m = 1/2. 
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FIG. 3. (a)Scaled plateau LA for ,h = 0.4, (b)central 
charge c and critical exponent rj, plotted versus the NNN 
coupling J2. 



IV. MAGNETIZATION PLATEAU 



We consider the magnetization plateau at m = 1/2. 
The plateau length A = E(M + 1) + E(M - 1) - 2E(M) 
is one of useful order parameters || to investigate the 
boundary between the gapless and plateau phases. Since 
A is the low-lying energy gap, it should obey the relation 
A ~ l/L in the gapless phase. The scaled plateau LA for 
several L is plotted versus J2 with fixed J\ to 0.4 in Fig. 
3 (a). It suggested that a gapless-gapful transition oc- 
curs at J2 ~ 0.2. To clarify the feature of the transition, 
we investigate the central charge c of the conformal field 
theory (CFT) Q and the critical exponent r\. r\ is de- 
fined by the asymptotic behavior of the spin correlation 
function (S£S~) ~ (-l) r r-''. CFT enables us to esti- 
mate c and 77 from the low-lying energy spectra of finite 
clusters, using the forms E(M)/L ~ e(m) — ttcv s /6L 2 
and A ~ Trv s r]/L (L — > 00), where v s is the sound veloc- 
ity which is the gradient of the dispersion curve at the 
origin. After some extrapolation to the infinite length 
limit, we show the results of c and 77 for J\ = 0.4 in 
Fig. 3(b). It justifies that the phase boundary is of the 
Kosterlitz-Thouless (KT) transition Jk}] with c = 1 in 
the gapless phase and rj = 1 at the critical point. Thus 
we determine the phase boundary as points with rj = 1 
in the J\-Ji plane. The result of the KT line is shown 
as solid symbols in Fig. 4 together with the first-order 
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boundary indicated as open symbols. Fig. 4 is a com- 
plete phase diagram at m = 1/2. The plateau phase is 
surrounded by the KT line and first-order line. The in- 
tersection of the two lines is expected to be a tri-critical 
point. The present analysis suggested that the plateau 
appears only in the rung-dimer phase. The rung-triplet 
phase reasonably has no plateau, because it is equivalent 
to the uniform 5 = 1 chain. 
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FIG. 5. Magnetization curves for (Ji, J2)=(0.5,0), 

(0.5,0.3) and (0.5,0.4). 

V. SUMMARY 

The antiferromagnetic spin ladder with NNN coupling 
is investigated by the exact diagonalization of finite clus- 
ters. It indicated the existence of the two magnetic 
phases; the rung-dimer and rung-triplet phases, not only 
for m = but also in the magnetic state. It is also 
found that the magnetization plateau possibly appears 
at m = 1/2 only in the rung-dimer phase. 



FIG. 4. Phase diagram at m = 1/2 including the plateau 
phase. The plateau appears only in the rung-dimer phase. 



A necessary condition of the presence of the plateau in 
general ID systems was rigorously given || by Q(S — m). 
Q is the periodicity of the ground state and S is the total 
spin of the unit cell. The present case must hold Q = 2 
in the plateau phase at m = 1/2. It suggests that the 
frustration stabilizes the structure where the singlet and 
triplet rung bonds are alternating, as is in the case of the 
zigzag ladder. @|| 

Finally we present the magnetization curves for 
(Ji,J 2 )= (0.5,0), (0.5,0.3) and (0.5,0.4) in Fig. 5. They 
were obtained by the size scaling in Ref. j|] applied to the 
calculated energy spectra of finite systems up to L = 16. 
The plateau clearly appears at m = 1/2 in the latter two 
cases. 
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